Supercurrents in an atom-molecule gas in an optical ring lattice 
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Atom and molecule currents in a Fermi gas in the neighborhood of a Feshbach resonance are 
studied in a one-dimensional optical ring lattice by directly diagonalizing small models. A rotational 
analogy of flux quantization is used to show that fraction of the current is carried by particles with 
twice the mass of an atom, which suggests pairing and superfluidity. 
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Artificial condensed-matter materials constructed by 
loading low-temperature atoms in an optical lattice hold 
promise as a means to explore major outstanding ques- 
tions such as the origin of high-T c superconductivity, to 
realize what heretofore have been fundamental thought 
experiments in condensed matter physics, and maybe to 
carry out quantum computations [J The degree of 
control exerted by the experimenter in an optical lat- 
tice stands in stark contrast to condensed-matter sys- 
tems supplied by Nature. For instance, the interactions 
between the atoms may be controlled by varying a mag- 
netic field around a Feshbach resonance 0, 0, Q . Pho- 
toassociation [f| 0, [1] could be similarly employed [§, [l(| • 

In fact, though, Feshbach resonance and photoassoci- 
ation entail conversion of atom pairs into molecules, and 
vice versa [ill ]. In this paper we study flows in an in- 
terconverting atom-molecule system in an optical lattice 
by direct numerical solutions of small systems. In a close 
analogy with past investigations of flux quantization [l2| , 
we analyze a one-dimensional optical lattice bent to a 
ring and rotated. We find indications for one facet of su- 
perfluidity with fermions, namely, that currents are car- 
ried by pairs of atoms reminiscent of Cooper pairs. Com- 
plete pairing may be approached at the Feshbach reso- 
nance when the atom-molecule coupling is strong, and 
when the magnetic field is varied across the resonance the 
carriers of the current smoothly change between atoms 
and molecules. 

What constitutes superfluidity is an interesting ques- 
tion in its own right. In the experiments on superflu- 
idity in optical lattices [H, [l4[ the criterion has been 
that, once released from the lattice, the atoms from dif- 
ferent sites interfere. This, however, is an observation 
of off-diagonal long-range order, not necessarily of super- 
fluidity. A persistent current is not a failsafe indicator 
of superfluidity either. In a mesoscopic ring that holds 
electrons and is threaded by a magnetic flux, the quan- 
tum mechanical ground state may carry a finite current. 
From this viewpoint, experimental observations of persis- 
tent currents in mesoscopic rings (Tol l are, perhaps, not 
surprising. Instead, a frequently used criterion for su- 
perfluidity in fermion systems is anomalous flux quan- 
tization. Basically, in phenomena associated with the 
presence of the magnetic field, the current appears to be 
carried by particles with twice the charge of an individual 
carrier [l3, [H, [l^ • 



Given that direct numerical solutions are feasible only 
in small systems, studies of long-range off-diagonal order 
that would model the experimental procedures [HI, [3| 
are presently beyond our capabilities. We therefore look 
for anomalous flux quantization. The remaining prob- 
lem with neutral atoms is that there is no direct cou- 
pling between magnetic field and the center-of-mass mo- 
tion. In principle one could use laser-induced electromag- 
netic c oupling s to produce an effective gauge field on the 
atoms [18j, |l9|, [20( , but here we envisage using rotation of 
the atomic sample as a substitute for magnetic field. 

Suppose that the prospective superfluid is tightly con- 
fined to a torus with the circumference L, and that in 
the direction of the torus, rotated periodic poten- 
tial V(x — vt) is imposed, where v is the linear ve- 
locity along the torus. When one transforms the one- 
particle Schrodinger equations for the wave function 
ip{x,t) to a co-moving frame, £ = x — vt, r = t, and de- 
fines ip(x,t) = e ™(^^/2+<)/fi^(^ T ) ; in tnc transformed 
Schrodinger equation the only remnant of the motion of 
the potential V(x) is the twisting boundary condition 
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*(^,r), where the rotation phase is 



$ = mvL/h. The value of $ only matters modulo 2ir. 

Next assume that the optical-lattice potential is deep 
enough to warrant a tight-binding approximation. The 
Wannier function for each site k — 1,...,N S , may be 
introduced as usual, u(x — Xk), but the annihilation op- 
erator of, say, a boson at site k, bk, is chosen in such a way 
that the corresponding wave function is e~ tktp u(x — Xk), 
ip = &/N s being the phase change from one lattice site to 
the next. The purpose of our convention is to revert the 
twisting boundary condition to the usual periodic bound- 
ary condition. The cost is altered phases of the hopping 
matrix elements, from a real t to te^ 1 ^ , in the ensuing 
Hubbard type model. 

The same approach works even if multiple interacting 
and interconverting species are present. We focus on a 
model in which two types of fcrmionic atoms (annihila- 
tion operators Ck\ and c&j,) may combine on-site to a di- 
atomic bosonic molecule (&&)• The Hamiltonian becomes 
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- tM^ie^bl^bk +h.c). 



2 



+ S^2blb k +g^2(blc kl c kr +h.c). (1) 

k k 

Here g is the coupling strength for atom-molecule conver- 
sion and 5 is the detuning from the Feshbach resonance 
controlled by the magnetic field [21]. We have thereby 
formulated what is known as a two-channel theory, with 
atoms and molecules as separate albeit coupled degrees 
of freedom [ll[ . We take the masses of the f and J, atoms 
to be equal, whereupon the mass of the molecules is twice 
the mass of the atoms and the two hopping phases satisfy 
ifM = Z^A- Different-mass atoms lead to incommensu- 
rate situations that do not have magnetic-field counter- 
parts, but we do not go there in this paper. Note that 
Hamiltonian (JTJ is the same as that for a ring lattice 
threaded by a magnetic field B = 2mvc/qL (22J. 

In the present model the conserved particle number 
N = J^k is a sum over the lattice sites of the local 

particle numbers N k = 4t c fct + 4i c *4 + 2b\bk- N k thus 
is a locally conserved quantity whose value only changes 
as a result of transport. This gives the identification of 
the operator for the current from the site k to the site 
k + 1 as 

4-fc+i = **A(e lVA 4 +1)tr c fe)(7 -e- lVA cl a c k+1 . a ) 

<r=T,4. 

+ 2it M (e iVM bl +1 b k -e- i ^blb k+1 ). (2) 

If the system is in an energy eigenstate in which there is a 
steady flow of atoms and/or molecules along the lattice, 
the value of the current I = (Ik^k+i) and the energy of 
the state E($ A , $m) regarded as a function of the phase 
parameters satisfies 

1 / dE n dE \ 

Since E is a period function of <& with period 2ir, so does 
/ as shown by Eq. ^ [23|. Here we may also compute 
the current numerically from Eq. |2]) and easily separate 
it to atomic and molecular components, but Eq. ^ is 
convenient for qualitative discussions. 

A way to compare the results for different system 
sizes is needed in order to ascertain that the outcome 
is not determined by limitations of the numerics. To this 
end we first note that for a nonrotating lattice, v = 
and ifA,M = 0, the ground state is obviously transla- 
tionally invariant and possesses a left-right symmetry: 

(4+l,T Cfe T> = (4'+i,T c fc'T> and (4+i,T Cfc t) = (4-i,T Cfe T> 
for all k and k' . The expectation values (4+i \°ki) are 
therefore real at v = 0. The motion of the lattice enters 
only through the phase parameters <pa,M- Assuming that 
perturbation theory in ifA,M is valid, Eq. @ shows that 
in the limit v ~ or tpA,M — the current is a linear func- 
tion of <pa,m — $A,m/N s . On the other hand, a sensible 
thermodynamic limit should be taken so that the num- 
ber of lattice sites N s and the invariant particle number 
N both tend to infinity in such a way that the average 
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FIG. 1: (Color online) Ground state energy E as a function of 
the rotation phase <3>a for g = (solid line), g/tA = 5 (dash- 
dotted line), and g/tA = 100 (dashed line). The energies 
corresponding to g/tA ~ 5 and g/tA = 100 have been scaled 
linearly so that the ranges of E coincide for all three g. The 
fixed parameters are 5 = 0, N 3 = 4, N = 2. 



occupation number per site N/N s remains a constant. 
Ground-state expectation values such as (4+i | c fct) m &y 
be viewed as some functions h^(N s , N/N s ), and in the 
limit N s — > oo with N/N s = constant they evidently 
tend to constants that depend only on the fixed ratio 
N/N s . Putting these considerations together, we see 
that, for given total phases &a and <I>m, the only possi- 
ble v-independent scaling of the current that leads to a 
nonzero and finite current in the thermodynamic limit is 
It = NJ. 

In this paper direct diagonalization based on the Lanc- 
zos algorithm [24[ is employed to find the eigenvalue (en- 
ergy) and eigenvector of the ground state of the Hamil- 
tonian (fTJ), and the results are then used to calculate var- 
ious observables at zero temperature. In this approach 
arbitrary parameter values may be accommodated with- 
out approximations, although only for small numbers of 
atoms and lattice sites. From now on we use the value of 
the atomic tunneling matrix element tA as the frequency 
scale. The number of f and J. fermions is always the same. 
Until further notice we also assume that the molecules do 
not tunnel, whereupon tui — and the parameters ipM 
and $m do not enter at all. 

We first study the resonance, S — 0, varying the cou- 
pling strength g. The dependence of E($a) on the phase 
&A is the weaker the larger is the coupling g. The current 
/ therefore decreases with g. With this in mind, we plot 
in Fig. [TJ the energies E as a function of $a for the cou- 
pling strengths g/t A — 0,5, and 100, scaling the energy 
linearly in such a way that the minima and the max- 
ima for each g coincide. The invariant particle number 
is N = 2, and the number of lattice sites is N s = 4. 

The main observation is that with increasing g the en- 
ergy changes from being a 27T periodic function of <fr A to 
having the period tt. The current, basically the derivative 
dE/d^A, behaves as if it were a sum of two components 
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FIG. 2: (Color online) Unpaired fraction U(g), a measure of 
the fraction of the current that is due to non-paired fermions, 
as a function of the coupling strength g. The fixed parameters 
are 8 = 0, N s = 8, N = 2. 



corresponding to the carrier masses m and 2m, with the 
mass 2m taking over as the coupling strength is increased. 
In analogy with the earlier studies of the influence of the 
magnetic flux [l2|, [H, H3] , we interpret this to mean that 
some of the fermions are paired, suggesting superfluidity. 
Moreover, the fraction of fermions appearing as pairs in- 
creases with increasing coupling strength. 

We next quantify the fraction of the current that is due 
to unpaired fermions. One might think that Fourier anal- 
ysis of E(<&a) gives straightforward answers, but there 
are problems. For instance, even with g = and no pair- 
ing at all, there is still a Fourier component in E($>a) 
with the period tt. We therefore use the ratio of the 
Fourier component of the current with the period 27r to 
an estimate of the current, normalized so that the value 
at g = equals unity, as the measure of the unpaired 
fraction. Specifically, we first find 

f( a ) = sin$J($; g )| 

ny> max<f[J($;j)] - min$[J($; g)] ' {> 

and define U(g) = f(g)/f{0). A representative unpaired 
fraction U is plotted in Fig. [2] as a function of the cou- 
pling strength g, for 8 = 0, N = 2, N s = 8. Unpaired 
fractions larger than unity are an unfortunate artifact of 
our definition, but the qualitative result is clear: With 
an increasing coupling strength the fraction of unpaired 
current carriers tends to zero. 

To investigate the effects of the necessarily small size 
of the system on the results we draw in Fig.[3]the current 
scaled by the number of lattice sites, It(&a) — N s I($a)> 
as a function of the rotation phase &a for the system 
sizes N s = 4 and N s = 8, both with N/N s = 1/2. We set 
8 = 0, and choose gjtA — 100. If It were an invariant in 
the thermodynamic limit, in the limit of large N s such 
curves should overlap. They do so quite well already in 
these small systems. Our interpretation is that, suitably 
scaled, the results for N s — 8 and N = 4 already rep- 



FIG. 3: (Color online) Scaled current as a function of the 
rotation phase for lattices with N a = 4, N = 2 (solid line) 
and N s — 8, N = 4 (dash-dotted line). The actual currents 
have been multiplied by the factors indicated in the legend. 
This figure is for Feshbach resonance, 5 = 0, and the atom- 
molecule coupling g/t a = 100 ensures that the fermions are 
almost completely paired. 



resent well the physics of the limit of a large number of 
lattice sites. 

We have also looked into the effect of varying the atom 
number N for a fixed number of sites N s on the fermion 
pairing behavior. In our numerics there is not much room 
for variation in N, and in a small system even a single 
unpaired fermion may make a major difference. With 
these caveats, in the strong-coupling case j > t/i we 
have not detected a qualitatively significant effect of atom 
number on the fraction of paired fermions as long as the 
number of fermions is even. 

In our final example we turn to the variation of the 
current with the detuning 8, the situation of BEC-BCS 
crossover where the system at near-zero temperature 
equilibrium is expected to change from an ideal Fermi gas 
to a BCS superfluid and then to a BEC of the diatomic 
molecules when the detuning is scanned across the res- 
onance starting from a large positive value. In such a 
process the increasing fraction of molecules, the "bare" 
molecules denoted by bk in our theory, has been observed 
experimentally [25] in a nonlattice gas, and successful 
comparisons of two-channel mean-field theories with the 
experiments have been reported [111 [27], [28| . Our model 
automatically leads to the period tt for the molecular cur- 
rents and as such says nothing about the superfluidity of 
the molecules, so that molecular currents are somewhat 
off the theme of the paper. Nonetheless, we shall demon- 
strate the change in the nature of the current carriers 
from atoms to molecules in a lattice upon the crossover. 

In order to allow for molecular currents in the first 
place we need to choose a nonzero hopping matrix ele- 
ment tM, thus we set tM — tA, and pick g/t a = 4.2 so 
that about half of the current is expected to be paired 
according to Fig. [5] We plot in Fig. 2] both the total cur- 
rent and the atomic and molecular currents separately as 
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a function of the varying detuning 5. The system size is 
N s = 8, N = 4. The expected smooth transition between 
atomic and molecular currents is observed. 
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FIG. 4: (Color online) Total current (solid line), atomic cur- 
rent (dash-dotted line) and molecular current (dashed line) 
as a function of detuning in a lattice with N s = 8, N = 4, 
t\i/tA = 1, and g/tA = 4.2. 

Although our study is meant to be an in-principle theo- 
retical discussion, a few comments about experiments are 
in order. Toroidal traps have been realized already [2^ |. 
and there is a proposal to make an optical ring lat- 
tice As usual, the lattice potential results from in- 
terference between two light beams. Also as usual, the 
lattice potential may be made to move simply by taking 
the two beams to have different frequencies, here possibly 
by means of an optical delay line incorporating a (slowly) 
moving mirror. Analysis of the velocity distribution of 
the atoms has been discussed in Ref. [2(j. Next, the cou- 
pling strength for the usual 834 G Feshbach resonance in 



6 Li is large compared to the widths of the energy bands 
in any conceivable optical lattice. This implies that a 
single-band model such as ([TJ) is inadequate [3(| ■ A faith- 
ful realization of the Hamiltonian ([1]) may have to await 
further developments in photoassociation (3lj . which in 
principle permits a full control of the coupling strength. 
As the tunneling matrix elements are exponentially small 
in the mass of the particle, the case with tji 3> tM is 
probably the experimental default. However, tunneling 
matrix elements can in principle be controlled with light 
by using optical transitions [321 ] . Overall, numerous and 
severe challenges must be overcome before our schemes 
are amenable to experiments, but the obstacles seem to 
be a matter of technology. 

In sum, we have studied currents in an interconvert- 
ing atom-molecule gas in an optical ring lattice by direct 
numerical solutions of small two-channel models. The 
center-of-mass motion of a neutral atom or molecule does 
not directly couple to the magnetic field, so we envisage 
rotation of the ring as a substitute to magnetic field for 
realizing an analog of flux quantization in a supercon- 
ductor. We find flux quantization as if part of the cur- 
rent were carried by atom pairs, which constitutes indi- 
rect evidence for superfluidity. With an increasing atom- 
molecule coupling, an increasing fraction of the current 
is carried by atoms pairs. In analogy with the BEC-BCS 
crossover, variation of the magnetic field in the neighbor- 
hood of the Feshbach resonance leads to a smooth switch- 
ing of the current carriers from atoms to molecules. The 
experimental challenges with our schemes are severe, but 
appear to be purely technical. 
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